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Abstract. We study one of the most X-ray luminous cluster of galaxies in the REFLEX survey, RXC J 1504. 1-0248 (hereafter 
^SJ . R1504; Zd = 0.2153), using XMM-Newton X-ray imaging spectroscopy, VLT/VIMOS optical spectroscopy, and WFl optical 

imaging. The mass distributions were determined using both the so-called hydrostatic method with X-ray imaging spectroscopy 
^ and the dynamical method with optical spectroscopy, respectively, which yield Mjjj^- = (5.81 ± 0.49) x 10''* Mq and MjjJ™'"^ = 

^w* ■ (4. 17 ± 0.42) X IO'^'Mq. According to recent calibrations, the richness-derived mass estimates closely agree with the hydrostatic 

' and dynamical mass estimates. The line-of-sight velocities of spectroscopic members reveal a group of galaxies with high 

velocities (>1000kms"') at a projected distance of about rjjjjj- = (1.18 ± 0.03) Mpc south-east of the cluster centroid, which 
is also indicated in the X-ray two-dimensional (2-D) temperature, density, entropy, and pressure maps. The dynamical mass 
estimate is 80% of the hydrostatic mass estimate at ''J'qq". It can be partially explained by the ~20% scatter in the 2-D pressure 
O ■ map that can be propagated into the hydrostatic mass estimate. The uncertainty in the dynamical mass estimate caused by 

• the substructure of the high velocity group is ~14%. The dynamical mass estimate using blue members is 1.23 times that 

^ I obtained using red members. The global properties of R1504 obey the observed scaling relations of nearby clusters, although 

its stellar-mass fraction is rather low. 

^SJ , Key words. Cosmology: observations — Galaxies: clusters: individual: RXCJ 1504.1-0248 — Galaxies: clusters: intracluster 

^ ■ medium — Methods: data analysis — X-rays: galaxies: clusters — Galaxy: kinematics and dynamics 

' 1 . Introduction ments the X-ray analysis in probing the dynamical structure of 

a cluster and the consequences of this for the mass measure- 
ment. 



■ Simulations show that the formation of galaxy clusters is not a 



purely gravitational process. For instance, the galaxy velocity 

dispersion in clusters provides evidence of heating when com- Comparisons between X-ray and weak-lensing mass esti- 



- '—I 
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pared to the cold dark matter velocity dispersion normalized n^^tes suggest that there is a 10-20% deviations from hydro- 
to the results from the Wilkinson Microwave Anisotropy Probe static equilibrium at rjoo for relaxed clusters of galaxies, and 



and large-scale structure distributions (e.g. Evrard et al. 2008). yield controversial findings for disturbed clusters based on dif- 
5^ . Cluster mergers do not only disturb the X-ray appearance of fgrent observations (e.g. Zhang et al. 2008, 2010; Mahdavi 
the hot intracluster medium (ICM; e.g. Ricker & Sarazin 2001 ; gt al. 2008). The comparison between the X-ray hydrostatic 
Poole et al. 2006; Zhang et al. 2009; Bohringer et al. 2010), and optical dynamical mass measurements complements the 
but also affect the properties of the cluster galaxies (e.g. Sun et x-ray versus weak-lensing mass comparison, and can thus 
al. 2007). The ICM and galaxies react on different timescales ^.^.^^^ ghegk deviations from hydrostatic equilibrium. Lensing is 
during merging (e.g. Roettiger et al. 1999). Clusters with sim- strongly affected by projection effects because it is sensitive to 
ilar X-ray global properties can contain rather different galaxy all mass along the line-of-sight. Line-of-sight velocities of clus- 
populations (e.g. Popesso et al. 2007). There is a higher frac- ^r galaxies via the redshifts of cluster galaxies can probe sub- 
tion of blue galaxies in rich galaxy clusters at z > 0.2 than in structures along the line-of-sight. Even for undisturbed clus- 
local clusters (Butcher & Oemler 1978, the so-called Butcher- ters, the dynamically relaxed central regions are surrounded by 
Oemler effect). This effect may be due to episodes of star for- infalling regions in which galaxies are bound to the clusters but 
mation in a subset of cluster members driven by merging (e.g. equilibrium (e.g. Biviano et al. 2006; Rines & Diaferio 
Saintonge et al. 2008). The kinematic analysis thus comple- 2OO6). Optical spectroscopy can constrain the dynamical prop- 
erties out to large radii (e.g. Braglia et al. 2007, Ziparo et al. 
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lites. On the other hand. X-ray data can be used to derive two- 
dimensional (2-D) cluster morphology projected on the sky. We 
use these two tools together with the optical colour information 
to constrain substractures, and model causes of deviations from 
hydrostatic equilibrium. 

RXC J1504.1-0248 (hereafter R1504; Zd = 0.2153), a 
prominent cool-core cluster, is one of the most relaxed clus- 
ters in the REFLEX galaxy cluster survey (Bohringer et al. 
2001). Bohringer et al. (2005) investigated the central region of 
R1504 with Chandra, and found that the core region of R1504 
appears to be quite relaxed. Using optical data obtained from 
integral field spectroscopy and X-ray data from XMM-Newton, 
Ogrean et al. (2010) found that feedback from active galac- 
tic nuclei (AGN) can slow down the cooling flow to the ob- 
served mass deposition rate in R1504, if the black-hole accre- 
tion rate is of the order of O.5M0 yr"' with 10% energy out- 
put efficiency. The XMM-Newton data in the R1504 field cover 
about a two times larger radial range than the Chandra data, 
and the large effective area of XMM-Newton ensures high qual- 
ity photon statistics for a 2-D spectral analysis. Follow-up ob- 
servations of R1504 were also performed by the ESO Visible 
Multi-Object Spectrograph at the 8.2 m UTS of the Very Large 
Telescope (VLT/VIMOS) out to about Srgoo. In addition, we 
have at our disposal optical imaging data from the ESO Wide- 
Field Imager (WFI) in three bands and Subaru in two bands. 

We present our analysis of both the radial and 2-D gas 
and galaxy properties of R1504 based on XMM-Newton, 
VLTA'IMOS, and WFI data as follows. We describe the X- 
ray analysis and results in Sect.|2] and optical spectroscopic 
and imaging analyses and results in Sect.|3]and Sect.|4] We dis- 
cuss the systematic uncertainties in and interpretations of the 
results in Sect.|5] In Sect.|6] we summarize our conclusions. 
Throughout the paper, we assume that On, = 0.3, Qa = 0.7, 
and //() = 70kms"' Mpc"'. Thus, an extent of 1' at the dis- 
tance of R1504 corresponds to 209.8 kpc. Confidence intervals 
correspond to the 68% confidence level. 

2. X-ray data analysis and results 

2.1. XMM-Newton observations and data preparation 

The XMM-Newton observations (ID: 0401040101) were taken 
in 2007 with a thin filter in full frame (FF) mode for MOS 
and extended full frame (EFF) mode for pn. The fraction of 
out-of-time (OOT) eff'ects amounts to 2.32%, which we used to 
normalize the pn OOT product before subtracting it from the 
pn normal product. 

We applied an iterative 2cr clipping procedure (e.g. 
Sect. 2.1 of Zhang et al. 2006) to filter flares using the light 
curves in both the soft band (0.3-lOkeV) binned in 10 s inter- 
vals and the hard band (10-12keV for MOS and 12-14keV 
for pn) binned in 100 s intervals. There are 34.9 ks clean data 
for MOS 1,34.5 ks for MOS2, and 3 1 .6 ks for pn left. We gener- 
ated a list of bright point sources detected with the XMMSAS task 
edetect_chain applied to five energy bands: 0.3-0.5 keV, 
0.5-2 keV, 2-4.5 keV, 4.5-7.5 keV, and 7.5-12 keV. All con- 
firmed sources were subtracted from the event lists using a 
radius of 25", which encloses ~80% of the flux of the point 



sources according to the XMM-Newton point spread func- 
tion (PSF). A weight column was created with the XMMSAS 
evigweight command, which accounts for the vignetting cor- 
rection for off'-axis observations. 



2.2. Temperature and surface brightness distributions 

Following Sect. 2.3 of Zhang et al. (2010), we measured 
the X-ray flux-weighted centroid to be RA=15*'04'"07':801, 
5=-02°48'ia.'29 (J2000). We applied the double-background 
subtraction method developed for clusters at similar redshifts 
by Zhang et al. (2006, the so-called DBS II method therein). 
We chose the XMM-Newton blank sky accumulations in the 
Chandra Deep Field South (CDFS), which used the medium 
filter and were also acquired in FF/EFF mode for MOS/pn as 
R1504, as background observations. The count rate in the hard 
band was used to normalize the CDFS observations before sub- 
tracting them from the target observations. The cluster X-ray 
emission was confined to be within R < 6', so that we use the 
annulus region outside R -W to probe the residual background 
and subtract it, taking into account the area coiTection for gaps 
and point sources. We chose the bin size of each annulus for 
the spectral extraction to be >0'.5 and to have about 4000 net 
counts in the 2-7keV band. The former ensures a redistribution 
fraction of less than 20% of the flux obtained within each an- 
nulus, and the latter ensures no more than 20% errors in the 
temperature measurements. Fig.[T]shows the X-ray spectra ex- 
tracted from the central <0'.5 circle. The spectra were fitted with 
the XSPEC software using the combined mekal*wabs model, 
in which the former describes the ICM emission, and the lat- 
ter describes the Galactic absorption with a Galactic hydrogen 
column density of n^. We freeze the cluster redshift to 0.2153 
and hh to the value of 6.08 x 10^° cm"- given by the LAB sur- 
vey (Hartmann & Burton 1997; Amal et al. 2000; Bajaja et al. 
2005; Kalberla et al. 2005) at the X-ray flux-weighted centroid. 

The left panel of Fig.|2] shows that the radial temperature 
distribution follows the universal profile of relaxed clusters 
(e.g. Markevitch et al. 1998; Vikhlinin et al. 2006; Zhang et 
al. 2006; Pratt et al. 2007). The observed strong cool core is 
also shown in Chandra data by Bohringer et al. (2005). The 
amplitude of the temperature distribution is, however, 20% 
lower than that of the Chandra measurements. Snowden et 
al. (2008) found that temperature measurements with Chandra 
are overestimated for those that are hotter than 5 keV owning 
to a calibration problem, which was corrected only recently 
with Chandra CALDBv3.5.2. The temperature overestimation 
made by analyses of previous Chandra CALDB thus accounts for 
the amplitude difference between the temperature distributions 
measured with Chandra by Bohringer et al. (2005) and with 
XMM-Newton in this work. The right panel of Fig.|2]shows the 
surface brightness distribution. We carried out the PSF decon- 
volution and deprojection as detailed in Sect. 2.3 of Zhang et al. 
(2007), in which the XMM-Newton PSF matrices were created 
according to Ghizzardi (2001). 
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2.3. Mass distribution using tlie liydrostatic method 

We assumed that the ICM is in hydrostatic equilibrium within 
the gravitational potential dominated by dark matter, and that 
the matter distribution is spherically symmetric. The mass 
could thus be calculated from the X-ray measured radial den- 
sity and temperature distributions of the ICM by 

1 d[n,{r)kT{r)] ^ GM(< r) 

in which ji - 0.62 is the mean molecular weight per hydro- 
gen atom, and k is the Boltzmann constant. The mass estimate 
was computed from a set of input parameters fi, ngOj, and 
(i - 1 , 2) representing the double-y6 electron number-density 
profile n,(r) = neOi(l + rV;^)-!-^/^ + n,02(l + r^r^T'-^P, and 
Pi (i - 1, ...,7) representing the deprojected temperature pro- 
file, T(r) = P3exp[-(r - PifP^^] + PeCl + r-p--)-''' + P^. 
We propagated the uncertainties in the electron number-density 
and temperature measurements using Monte Carlo simulations 
to compute the mass uncertainty. The cluster radius (e.g. 
rsoo) was defined as the radius within which the mass density is 
A (e.g. 500) times of the critical density at the cluster redshift, 
p^{z) = 2,Hl{%nG)-^E\z), in which E\z) = Qm(l +1? + Q.^ + 
(1-0„-Qa)(1+z)2. 

The cumulative mass distribution using the hydrostatic 
method is shown in Fig. [3] as filled squares. The cluster ra- 
dius is r^gy- - (1.18 ± 0.03) Mpc, and the total mass estimate 
is M^,E. ^ (5 82 + 0.49) X 10"*Mo as listed in Table[T] The 
mass profile can be well-fitted by an NFW model (Navarro et 
al. 1997) with a concentration parameter of C500 - 2.93 + 0. 14. 

2.4. Two-dimensional spectrally measured maps 

We adopted the adaptive-binning methods developed by 
Sanders (2006) and Cappellari & Copin (2003; see also Diehl 
& Statler 2006), respectively, to define the masks for the 2-D 
analysis. We followed the steps in Sect. 5 of Zhang et al. (2009), 
except for two changes: (i) We adopted a threshold of signal-to- 
noise ratio (S/N) of 30 instead of 45 to ensure an adequate num- 
ber of bins in the mask, (ii) We used the double-background 
subtraction strategy (i.e. DBS II) since there is a sufficient outer 
area within the XMM-Newton field of view (FoV) to model the 
residual background. 

The X-ray spectrally measured 2-D maps using the Sanders 
(2006) mask are shown in Fig.|4] Those using the Cappellari 
& Copin (2003) mask appear in Fig. lA.ll At a projected dis- 
tance of ~ r^jp south-east of the centroid, the surface bright- 
ness shows a clump with its value 2.25 times the azimuthal 
average of (2.4 ± 0.4) x 10""' counts s"' arcmin"^, at the same 
projected distance. This enables us to carry out the 2-D spectral 
analysis toward larger radii in the south-east direction than the 
other directions. In the temperature map, we also observe a hot 
stripe at ~ 0.3r^,jQ' west of the centroid. 

We followed Zhang et al. (2009) to calculate the scatter in 
the fluctuations in the 2-D maps and analyse the ICM substruc- 
ture. In Fig. |5] the temperature, electron number-density, en- 
tropy, and pressure distributions in the 2-D maps are shown as 
a function of projected distance from the X-ray flux-weighted 



centroid. Their non-parametiic locally weighted regression 
(e.g. Becker et al. 1988) was used as the mean profile. The 
differential scatter and error were calculated from the area- 
weighted absolute fluctuations with all bins at that projected 
distance. The temperature, electron number-density, entropy 
and pressure maps have scatters of smaller than ~15%, 30%, 
20%, and 25%, respectively. The scatter in the temperature 
measurements increases with radius, and is ~ 12% at the radii 
beyond 0.3r^Qg . In contrast, the electron number-density mea- 
surements appear to have a smaller scatter at large radii, which 
may reflect the different relaxation timescales of the ICM tem- 
perature and electron number density. 

The 2-D maps using the Cappellari & Copin (2003) mask 
demonstrate a mild elongation along the east-west major axis 
on the scale of ~ r™ • . 



3. Optical spectroscopic data analysis and results 

3.1. VLT/VIMOS observations and data preparation 

Multi-object spectroscopy in the R1504 field was performed 
with the VLT/VIMOS spectrograph using the low-resolution 
blue grism (LR-blue) with the OS-blue filter, which samples 
the [3700-6700] A wavelength range with A/AA = 180 (PID: 
077.A-0058). The observations were caiTied out in service 
mode, which include a bias, dark, and flat calibration sequence 
with the last of these using the same set-up as the observations. 

The cluster RI504 was observed with air masses of be- 
tween 1.124 and 1.172. The seeing was less than r.'5 for two 
pointed observations and > 2" for one pointed observation. 
Individual targets for the spectroscopic follow-up were selected 
based on their /-band magnitudes in the pre-imaging data. In 
total, 784 slits were placed in the 12 masks spread over three 
pointed observations. 

Data were reduced with the VIPGI pipeline (Scodeggio 
et al. 2005). Redshifts were computed by comparing the ob- 
served spectra to galaxy and stellar templates with the EZ 
software (Garilli et al. 2010) and improved with customized 
tools from Verdugo et al. (2008) by fitting a Gaussian pro- 
file to each of the observed strong spectral features. Not 
all lines, namely [On] ^1/13726,3729 AA, [Can K] i3934A, 
[Can H] /13968A, H6 /14102A, Hj3 /14861A, G-band (~ 
^4304 A), [Oni] ^1/14959,5007 AA, Mg lb (~ 5050 - 5430 A), 
and [Fe VI] /15335 A, are present in each spectrum. The red- 
shift determined from each spectrum is the mean of the indi- 
vidual shifts for all visible line features. The error is the stan- 
dard deviation in these shifts, which takes into account system- 
atic errors in the wavelength calibration and differences in the 
S/N. The average error in the individual redshift estimates is 
Az = 0.00076. 

The R1504 field overlaps with the stellar stream of the glob- 
ular cluster Palomar-5 (e.g. Odenkirchen et al. 2001). There are 
thus a large number of white- and red-dwarf stars contaminat- 
ing the spectroscopic sample. Interlopers were removed by ap- 
plying the member selection procedure in Sect. 13.21 
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3.2. Spectroscopic member selection 

In hierarchical structure-formation scenarios, the spherical in- 
falling model predicts a trumpet-shaped region in the diagram 
of line-of-sight velocity versus projected distance, the so-called 
caustic (e.g. Kaiser 1987). The boundary defines galaxies in- 
side the caustic as cluster members and those outside the caus- 
tic as fore- and background galaxies. 

Several methods can select cluster members according to 
their line-of-sight velocities. A constant gap or gaps weighted 
by projected distances can separate member galaxies from fore- 
and background galaxies (e.g. Beers et al. 1990; Girardi et al. 
1993; den Hartog & Katgert 1996; Katgert et al. 1996, 2004; 
Popesso et al. 2005a; Biviano et al. 2006). Alternatively, an 
adaptive kernel method, sensitive to the presence of peaks, has 
been widely used to select members (e.g. Pisani 1993, 1996; 
Fadda et al. 1996; Diaferio & Geller 1997; Diaferio 1999; 
Geller et al. 1999; Rines & Diaferio 2006). 

We modified the adaptive kernel method of Diaferio (1999) 
for medium-distance clusters (e.g. Harrison & Noonan 1979; 
Danese et al. 1980), and used it to identify the interlopers as 
follows. Galaxies with redshifts \cz - czcil < 4000 km s"' were 
preliminarily selected and located in the (r, v) diagram, where 
r is the galaxy transverse separation from the cluster centroid, 
and V is the line-of-sight velocity of the galaxy relative to the 
cluster At the redshift of R1504, we modified the calculation 
of (r, v) to r = Da sin and v = {cz - czcicos0)/(l + Zd), 
where D^, 9, z, and Zd are the angular diameter distance of the 
cluster, angular separation of the galaxy from the cluster cen- 
troid, galaxy redshift, and cluster redshift. The angular diam- 
eter distance is defined as Da = c//q'(1 + Zci) ' E^\z)dz- 
We scaled r and v to ensure approximately equal weights in 
r and v, and computed the 2-D adaptive density distribution 
f(r, v). The boundary of the caustic is defined by f{r, v) - k, 
in which k is determined by minimizing \{vI^^}k,r - 4(v')r|~. 
Here {v^sc}k,r - Jq Jr(r)ip{r)drl J^^ tp{r)dr is the escape ve- 
locity, R and (v-)r are the mean projected distance and velocity 
dispersion of the members, and ip(r) = J f{r, v)dv. In addition, 
^(r) - min{|vu|, |vd|), which represents the minimum of the up- 
per and lower solutions of f{r, v) - k. 

Fig.|7] shows the histogram of the line-of-sight velocity 
and the diagram of line-of-sight velocity versus projected dis- 
tance of galaxies. The 53 galaxies inside the caustic boundary 
are considered as members. Our further dynamical analysis is 
based on these 53 spectroscopic member galaxies. 

3.3. Mass distribution using tine caustic method 

We applied the biweight estimator developed by Beers et al. 
(1990) to measure the redshift and velocity dispersion, in which 
the errors are estimated by means of 1000 bootstrap simula- 
tions. The cluster redshift was assumed to be the biweight es- 
timator of location of the recessional velocities, cz,. The ve- 
locities in the cluster reference frame Viest,/ were derived from 
the recessional velocities cz, using the formula Viest.i - {czi - 
cZci)(l + Zci) The velocity dispersion was given by the bi- 
weight estimator of scale of these Vrest,; values. With all 53 spec- 



troscopic members, the measured redshift is (0.2165 ± 0.0005), 
and the velocity dispersion is (1 132 + 94) km s"' . 

The escape velocities of galaxies extracted from the ampli- 
tude of the caustic scale with the gravitational potential of the 
dark matter halo (e.g. Diaferio & Geller 1997; Diaferio 1999; 
Rines & Diaferio 2006) 

1 f 

GM(< r) = - Jf(x)dx. (2) 

2 Jo 

The error in the cumulative mass estimate within r,- is defined 
as SMi = Tj'j=\ |2'«>5^('"y)/^('"y)l, in which mj is the mass 
estimate in the shell, 6^{r)/^{r) = a:/ max{/(r, v)), and 
max{/(r, v)) is the maximum value of f(r, v) at fixed r. Since 
the caustic of a spherical system is symmetric, the minimum of 
the upper and lower caustics is taken as the caustic amplitude 
at each radius, which excludes interlopers more effectively than 
taking the average of the upper and lower caustics as the caus- 
tic amplitude. The caustic method does not assume dynamical 
equilibrium, and is thus able to measure cluster masses out to 
large radii. The mass distribution using the caustic method is 
shown in Fig. |3] and implies that r^™'*'''= = (1.06 + 0.04) Mpc 
(Table©. 

3.4. Two-dimensional kinematic structure 

As clusters of galaxies formed at moderate look-back times, 
they often contain substructures caused by the accretion of and 
merging with smaller systems. The quantification of substruc- 
tures is thus important to assess the reliability of different mass 
measurements and infer the dynamical history of the cluster. 
It is, however, diflicult to quantify substructures at large radii 
with XMM-Newton data because of the high background level 
relative to the faint X-ray emission near /"^gy '. Optical spectro- 
scopic data provide a promising alternative means of quantify- 
ing substructures at all radii. 

Dressier & Shectman (1988, DS hereafter) developed a sta- 
tistical method for detecting substructures in galaxy clusters 
related to systematic deviations from the average spatial and/or 
velocity structure, which we summarize as follows. The devia- 
tion in the local velocity mean and dispersion from the global 
mean and dispersion of the cluster is defined as 



Pds = y~~2 [(^local - V')^ + (O-Iocal - 0")^] . (3) 

The local line-of-sight velocity and dispersion, i'locai and crjocai, 
are measured with its A^iocai - 10 nearest neighbours in pro- 
jection. The global line-of-sight velocity and dispersion, v - 
64950 km s"' and cr = 1132kms"', are given in Sect. 13.31 as 
measured for all (i.e. N^n = 53 for R1504) spectroscopic clus- 
ter members. A cumulative deviation Aqs = Tjf^i Pdsj is of the 
order of A^aii when the line-of-sight velocity obeys a Gaussian 
distribution and local variations are only random fluctuations, 
but differs from A^aii when the line-of-sight velocity distribu- 
tion deviates from a Gaussian. Therefore, the Ads statistic is 
sensitive to the presence of substructures. A higher Aqs value 
indicates a higher possibility of containing substructures. 

In Fig. 3] we display the 53 spectroscopic galaxies by cir- 
cles with their radii proportional to exp[ppg]. The blue, red, 
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and green colours denote cluster galaxies that have spectro- 
scopic follow-up data with their clustercentric velocities to- 
ward the observer greater than 1 000 km s " ' , away from the ob- 
server greater than 1000 km s"', and smaller than 1000 km s"'. 
The DS test yields Aos.obs = 86.8, and reveals a clustering of 
galaxies with large pus values at a projected distance of ~ r^jy • 
south-east of the cluster centroid. This may be indicative of a 
relatively high-velocity (> 1000 km s"') group associated with 
the main cluster. 

To estimate the robustness of the Ads statistic, we car- 
ried out hmc = 10,000 Monte Carlo simulations by randomly 
shuffling the line-of-sight velocities of the member galaxies. 
This procedure erases any true correlations between the line- 
of-sight velocities and positions. The probability of having sub- 
structures can be estimated by Pus = ^j/^MC, in which 
6j - Adsj when Aosj > Aos.obs and 6j = otherwise (e.g. 
Hou et al. 2012). Therefore, systems with significant substruc- 
tures show low-f DS values. The Monte Carlo realization of the 
Ads statistic yields f ds - 0.06 for R1504, which only cor- 
responds to a ~ 10% false detection according to Hou et al. 
(2012). This f DS-value thus represents a non-negligible possi- 
bility that the high velocity group exists. The histogram of the 
line-of-sight velocities of possible member galaxies of the high 
velocity group is shown as (red) dashed lines in the left panel 
ofFig.El 



4. Optical imaging data analysis and results 

4.1. WFI observations and data preparation 

The cluster R1504 was observed with the ESO/MPG 2.2 m 
telescope with the WFI (Baade et al. 1999) in the B/123-, 
V/89-, and /?<;/162-bands (hereafter B-, V-, and -bands) dur- 
ing three runs in June 2009, 2010, and 201 1 mostly under non- 
photometric conditions. The WFI has a 4 x 2 mosaic detector 
of 2k X 4k CCDs, which has a pixel size of a.'238 and a FoV 
of 34'x33'. We reduced the data and co-added the single expo- 
sures using the THELI pipeline (Erben et al. 2005). The coad- 
ded B-, V-, and /?-band images have total exposure times of 
9238 s, 8317 s, and 11035 s with 1'.'27, 1'.'21, and ^.'96 seeing. 
Fig.[8] shows the WFI B-band image. 

Since the observing conditions of the imaging varied 
slightly, photometric calibration of the WFI 7?-band data was 
tied to the shallower SDSS r-band data in the same field. The 
B- and V-band images were calibrated against the 7?-band by 
fitting the observed main sequence of stars with that of a photo- 
metrically calibrated field. Cross-checks were applied by com- 
paring predicted with observed colours for galaxies with spec- 
troscopic redshifts using procedures typical of photometric red- 
shift codes. The absolute calibration of the photometric zero 
point is more accurate than 0.1 magnitude. The internal colour 
calibration is more accurate than 0.04 magnitude. 

Photometry was carried out with SExtractor (Bertin & 
Arnouts 1996) in dual mode with the 7?-band image set as the 
detection frame. Galaxy colours were measured in an aperture 
of 2'.' 5 diameter. 



4.2. Photometric member selection 

The broad B- and /?-bands can probe the spectral range across 
the 4000 A break for galaxy clusters at z ~ 0.2 - 0.3, which can 
be used to characterize the bulk of stellar populations. Early- 
type galaxies display red colours and reside in a naiTow band, 
the so-called red sequence (e.g. Baum 1959), in the colour- 
magnitude diagram. The B-R colour versus /?-band magnitude 
is shown in the left panel of Fig.|9]for all galaxies within /"^qq , 
in which we highlight the 29 blue and 24 red spectroscopic 
members regardless of their clustercentric distances. 

We fit the red sequence with a robust biweight linear-fitting 
algorithm, which is negligibly aff'ected by skewness (e.g. Beers 
et al. 1990; Gladders et al. 1998). The best-fit relation is B-R = 
-0.0977? + 3.88 with a scatter of cr = 0.083 magnitude. 

We define red members as galaxies within 3cr of the best- 
fit red sequence. Galaxies that are more than 3cr bluer than the 
red sequence are considered as blue galaxies, whereas galaxies 
that are more than 3cr redder than the red sequence are dis- 
carded because they are most likely background galaxies. We 
also discard galaxies bluer by more than 1 .5 magnitude than the 
red sequence to reduce foreground contamination. These limits 
are marked in Fig.|9] 

The comparison between the galaxy number counts in the 
R1504 field and the Garching-Bonn Deep Survey (GaBoDS, 
Erben et al. 2005; Hildebrandt et al. 2006) fields shows that 
our photometric sample is 100% complete down to R - 
22 magnitude. This limit coiTesponds to Mr x -18 magnitude 
for galaxies of different colours atZci = 0.2153, and is indicated 
by an oblique curve in Fig.|9] The CLASS_STAR parameter pro- 
vided by SExtractor is completely reliable down to R - 
21.5 magnitude and 90% reliable down to R - 22 magnitude. 
We thus limit the analysis to the galaxies brighter than R - 
22 magnitude. The limiting magnitude is R - 24.44 magnitude 
at 5cr, which was calculated in an aperture of 1 " diameter. The 
stellar contamination should be quite low in the photometric 
sample because of the high quality of the WFI imaging. 

4.3. Galaxy luminosity function and stellar mass 

To determine the photometric properties of the cluster popula- 
tion, we corrected for contamination by both fore- and back- 
ground galaxies along the line-of-sight by applying a statis- 
tical subtraction using a subset (2.25 deg') of the GaBoDS, 
whose photometric data are complete in all bands. We applied 
the same magnitude and colour cuts to the GaBoDS fields as 
we did to the R1504 field, and constructed the galaxy number 
counts as a function of magnitude for red and blue galaxies, 
respectively. The galaxy number counts in the GaBoDS fields 
were normalized to the R1504 sky area and subtracted from the 
corresponding galaxy number counts in the R1504 field. The 
size of the magnitude bin used to construct the galaxy num- 
ber counts was chosen to be large enough to ensure that there 
were positive counts in all magnitude bins after subtracting the 
background. The errors are a combination of the Poisson sta- 
tistical error and cosmic variance. The latter was measured to 
be the standard deviation in the galaxy number counts across 
the GaBoDS fields. The Poisson statistical error dominates the 
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error budget at inner radii because there are relatively fewer 
objects, whereas cosmic variance dominates at large radii. 

The right panel of Fig.|9]shows the 7?-band galaxy luminos- 
ity functions for red (red), blue (blue), and all (black) members 
within r^QQ'- We derived the best-fit Schechter (1976) functions 
for red and all members, respectively, with a minimization 
algorithm. The red members show a shallower slope at the faint 
end than that of all members. As listed in Table[T] both the 
characteristic 7?-band magnitude, R*, and the slope at the faint 
end, a, agree with that of nearby clusters (e.g. Christlein & 
Zabludoff 2003; Popesso et al. 2005b; Durret et al. 2011). No 
fit was possible for blue galaxies. 

The brightest cluster galaxy (BCG) is very blue (B - R = 
0.98) and bright (R - 16 magnitude). Since the BCG cannot 
normally be described by the Schechter function, we explicitly 
removed the BCG from the photometric sample in the above 
calculation of the galaxy luminosity function. The total /?-band 
luminosity of the cluster was computed by explicitly adding 
the BCG luminosity to the integration of the galaxy luminosity 
function. We tested how the BCG colour affects the measured 
optical luminosity by adding the BCG luminosity to the red and 
blue populations, respectively. The resulting uncertainty in the 
7? -band luminosity is smaller than 3%. The measured ^-band 
optical luminosity within r^^- is fisted in Table[T] 

We adopted the stellar mass-to-optical light ratio from 
Table? of Bell et al. (2003), logiMJlg) ^ a + b(B - R), in 
which a - -0.523, b - 0.683, B - R is the rest-frame colour, 
and a "diet" Salpeter (1955) stellar initial mass function (Bell & 
de Jong 2001) is adopted. We show the cumulative stellar mass 
profiles for red and blue galaxies, respectively, in Fig.|3] The 
stellar mass estimates for blue, red, and all photometric mem- 
bers within r^yy are given in Table[T] The cluster R1504 has 
a rather low stellar-mass fraction among nearby clusters (e.g. 
Biviano & Salucci 2006), although there is an agreement given 
the large scatter for nearby clusters (e.g. Andreon 2010; Zhang 
et al. 2011b). The contribution to the total stellar mass budget 
associated with intra-cluster light (e.g. Zibetti et a. 2005) could 
reach 10-20% level at the mass scale of galaxy groups (e.g. 
Gonzalez et al. 2007). The difference in the /?-band luminos- 
ity caused by adding the BCG luminosity to the blue and red 
populations, respectively, is within 3% for R1504. Therefore, 
both the intra-cluster light and BCG play non-dominant roles 
to explain the low stellar-mass fraction. 

4.4. Blue-galaxy fraction 

Understanding rapid galaxy-morphology transformation is es- 
sential to model the mass assembly history of galaxy clus- 
ters. Red galaxies with low levels of ongoing star formation 
are assumed to be the descendants of blue star-forming galax- 
ies that have been accreted from the surrounding filamentary 
structure and quenched by specific processes operating within 
the cluster environment (e.g. Braglia et al. 2007; Verdugo et al. 
2012). How spirals were transformed into ellipticals remains 
unclear (e.g. Poggianti & Wu 2000). The efficiency of dif- 
ferent transformation processes varies with environment. Ram 
pressure stripping, for instance, is more eff'ective in the cluster 



core where there is dense ICM (e.g. Fujita & Nagashima 1999; 
Quilis et al. 2000). Galaxy-galaxy merging, on the other hand, 
is more effective in infalling groups at the cluster outskirts (e.g. 
Kauff^mann et al. 1999). 

In the left panel of Fig.[10l we present the number density 
profiles of the red and all member galaxies, respectively, as a 
function of projected distance from the X-ray flux-weighted 
centroid. Similar to the findings for other clusters (e.g. Lin et al. 
2004), the red population in R1504 is well-fitted by a projected 
NFW profile (Navarro, Frenk & White 1997; Bartelmann 
1996), in which the reduced is 1.03. The fit of a projected 
NFW profile to all members is relatively poor with = 2.24. 
The right panel of Fig.[TO] displays the blue-galaxy fraction. 
Within a projected distance of 700 kpc, where the ICM is dense, 
the galaxy population is dominated by red galaxies with a frac- 
tion larger than 70%. Beyond 700 kpc, the blue-galaxy fraction 
first increases rapidly, reaching 60% at ~ '"500 > ™d then stays 
almost constant in the outskirts. 

The large error bars in the blue-galaxy fraction at large 
radii take into account the uncertainty in the background. The 
colours of the field galaxies have their own variation. One can- 
not exclude that the high blue-galaxy fraction is coincidently 
due to a larger field-galaxy contamination. A reliable mea- 
sure of the blue-galaxy fraction free of background contami- 
nation requires extensive spectroscopy. Nevertheless, the high 
blue-galaxy fraction beyond r^^- may indicate a critical dis- 
tance where specific processes in clusters such as starvation 
(e.g. Balogh et al. 2000) start to affect the properties of in- 
falling galaxies. This distance may, however, depend on the 
environment (e.g. Urquhart et al. 2010). Urquhart et al. (2010, 
Figure 10 therein) pointed out an excess of extremely blue 
galaxies in massive groups, in which active starburst galaxies 
are driven by galaxy-galaxy interactions in the group environ- 
ment. The radial blue-galaxy fraction and presence of the high 
velocity group in R1504 support this interpretation. 

4.5. Richness-derived mass estimates 

There are a number of approaches for calculating the cluster 
total mass based on photometric measurements (e.g. Hansen et 
al. 2005, 2009; Johnston et al. 2007; Reyes et al. 2008). 

We can iteratively estimate the mass from the richness 
distribution according to existing scaling calibrations between 
richness and mass estimates at a fixed overdensity. The cluster 
total mass is measured within the radius where the mass den- 
sity is 200 times the critical density in Johnston et al. (2007) 
and Hansen et al. (2009), but 200 times the mean density of 
the Universe in Hansen et al. (2005) and Reyes et al. (2008). 
To be consistent with the mass definition, we use the calibra- 
tions between richness and weak-lensing mass measurements 
in Johnston et al. (2007), eq. (26), and Hansen et al. (2009), 
eq. (10), respectively, to calculate the richness-derived mass es- 
timates. In these calibrations, the richness was measured down 
to 0.4L* in the SDSS /-band, which is equivalent to (/* + I) 
magnitude. Since there were no /-band data for R1504, we used 
the /?-band derived richness down to (R* -1- 1) to approximate 
their /-band derived richness down to (/* + 1). Nevertheless, 
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the characteristic magnitude could be taken as a characteristic 
length-scale for galaxy formation at a similar cosmic epoch. 
Therefore, galaxies brighter than (/* + 1) overlap significantly 
with those brighter than {R* + 1). The richness-derived mass 
estimates (see Table[T]) agree well with the hydrostatic and dy- 
namical mass measurements. 

5. Discussion 

5.1. Uncertainties in the dynamical mass estimates 

To quantify the uncertainties in the dynamical mass estimate 
due to substructures and infalling blue galaxies, we need a 
method to derive the dynamical mass estimate based on a small 
number of member galaxies. A different method from the caus- 
tic one, which was detailed in Sect. 3 of Biviano et al. (2006) 
based on a sample of simulated clusters, serves this purpose, 
and is summarized as follows. 

We measured the velocity dispersion (cr^ p in Biviano et al. 
2006) with a number of member galaxies defined according to 
Sect. 15.1.11 and Sect. 15.1.2! respectively. An initial estimate of 
the mass was derived using eq. (2) in Biviano et al. (2006), 
namely My = A[crJ(10^ km s"')]^^ x 10'4/2-i£-i(z)Mo with 
A - 1.50 ± 0.02 under the assumption of cry - V3cr;,p. An 
estimate of the corresponding cluster radius, fy, was derived by 
following steps 7 and 9 in Biviano et al. (2006). Replacing the 
true quantities ry and cr.^ with their estimates ry and VScr^.p in 
Fig. 4 in Biviano et al. (2006), we obtained an improved es- 
timate of CTy. The final cluster mass estimate {M^^ hereafter) 
was calculated using eq. (2) in Biviano et al. (2006) with the 
improved estimate of cry. The mass error was derived by com- 
bining in quadrature the error in the velocity dispersion con- 
verted to mass, the additional error introduced by the uncer- 
tainty in Fig. 4 in Biviano et al. (2006), and the 10% mass sys- 
tematic according to the Biviano et al. (2006) sample. We com- 
puted MjQ°Q from MjyQ* with the NFW model, in which the con- 
centration parameter is given in step 7 of Biviano et al. (2006) 
as c^°^ = 4[o-a,p/(700kms-i)]-"-3°<'. 

5.1 .1 . Uncertainty due to substructures 

The high velocity group revealed by the DS test overlaps with 
the substructure in the 2-D X-ray maps. In the high-resolution 
weak-lensing shear map by Klein et al. (in prep.) based on 
Subaru observations, there is a main lensing peak at the BCG 
position additional to a secondary peak near the high velocity 
group. The elongation of the surface mass density follows the 
elongation shown in the X-ray 2-D maps. 

The dynamical mass measurement overestimates the total 
mass of a cluster with substructures (e.g. Biviano et al. 2006). 
We quantified the uncertainty in the dynamical mass estimate 
due to the high velocity group as follows. The velocity disper- 
sion is (1 132 + 94) km s"' based on all 53 spectroscopic mem- 
bers. Excluding the five galaxies belonging to the high velocity 
group, which are shown as large red circles at ~ r^j^- south- 
east of the cluster centroid in Fig.|4] the velocity dispersion is 
reduced to (1079 + 97) km s"' . As shown in Table[T]and Fig. [3] 
the total mass estimate is M^!^„ = (6.355 + 1.71 1) x IO'^^Mq 



based on the 53 members and Mf"?,„ = (5.452 ± 1.573) x 

500,48m ^ ^ 

IO'^^Mq based on the 48 members excluding the high velocity 
group, following the Biviano et al. (2006) method. The mass 
uncertainty caused by the substructure of the high velocity 
group is ~ 14%. Excluding the high velocity group improves the 
agreement between the X-ray hydrostatic and dynamical mass 
estimates. This highlights the importance of a proper analysis 
of the structures of evolving systems such as cluster of galaxies 
in order to obtain robust mass measurements. 

5.1 .2. Mass discrepancy using red and blue galaxies 

For a cluster in dynamical equilibrium, both blue and red 
galaxy populations behave similarly in the diagram of line- 
of-sight velocity versus projected distance, which reflects the 
same underlying mass distribution. For a cluster with a rich 
infalling population, the blue population tends to overestimate 
the total mass since the infalling members bias the velocity dis- 
persion toward high values (e.g. Carlberg et al. 1997). Fifty-five 
percent of the spectroscopic member galaxies are blue galax- 
ies. To quantify the uncertainty in the dynamical mass estimate 
due to infalling blue galaxies, we derived the velocity disper- 
sion and dynamical mass estimates for blue members and red 
members separately. The blue population has a velocity disper- 
sion of (1170 + 125)kms"', which yields a cluster mass es- 
timate of M^°„^y^,^ = (7.057 ± 2.373) x lO'^^Mg according to 
Biviano et al. (2006). The red population has a velocity disper- 
sion of (1096 + 156) km s"', which yields a cluster mass esti- 
mate of M™f^.^^ = (5.722 + 2.517) x lO'^'Mg. The dynamical 
mass estimate derived from red members, M^l^^.^^, is in closer 
agreement with the hydrostatic mass estimate than that using 
all members. The cluster mass estimate based on the blue pop- 
ulation, M^nnu, , is 1.23 times the mass estimate based on the 

' 500,blue' 

red population (Table[T]and Fig.|3]l. It indicates that R1504 may 
not contain a rich infalling population, which is also supported 
by the caustic comparison between R1504 and the simulated 
clusters (e.g. Figure 6 of Gill et al. 2005). 

5.2. Hydrostatic versus dynamical versus photometric 
mass estimates 

We derived the mass distribution using the hydrostatic method 
based on X-ray imaging spectroscopy and the dynamical 
method based on optical spectroscopy, independently (Fig.O. 
At ^"^00 ' mass estimate from the caustic method is (4.66 + 
0.47)x 1O''*M0, which is 80% of the X -ray hydrostatic mass es- 
timate (TablelTJ. The discrepancy between the hydrostatic and 
dynamical mass estimates is comparable to the mass discrep- 
ancy between the results obtained using red and blue galaxies. 
In addition, the 2-D pressure map shows a ~20% scatter, which 
can be propagated into the hydrostatic mass estimate. Despite 
the mass uncertainties, the X-ray hydrostatic mass estimate is 
slightly higher than the dynamical one at r^^g because the mass 
profile of the latter appears to be more concentrated than the 
former. 

As shown in Table[T] the richness-derived mass estimates 
according to recent calibrations (e.g. Johnston et al. 2007; 
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Hansen et al. 2009) agree well with the hydrostatic and dy- 
namical mass estimates. 

It is also important to compare the mass measurements at 
fixed overdensities (e.g. 500 and 200) as these masses are used 
in cosmological applications. At the overdensity of 500, the 
mass estimate from the caustic method is 72% of the hydro- 
static mass estimate. This is due to the less-concentrated mass 
distribution determined using the hydrostatic method than that 
determined from the dynamical method, which can be caused 
by the disturbed ICM in the cluster core because of AGN 
feedback as discussed in Sect. 15.41 At the overdensity of 200, 
the richness-derived mass estimates according to Hansen et 
al. (2009) and Johnston et al. (2007) are about 0.8-1.2 times 
the hydrostatic mass estimate and 1.0-1.4 times the dynamical 
mass estimate. 

5.3. Agreement with the scaling relations 

The X-ray properties of R1504 are tabulated in Table[T] The 
global temperature and metallicity are measured with the 
XMM-Newton spectra extracted from the (0.2 - 0.5)r^QQ- an- 
nulus. The X-ray luminosity was derived by integrating the 
surface brightness out to r^^\ We show the X-ray luminos- 
ity of the 0.1-2.4keV band and 0.01-100 keV (bol hereafter) 
band, respectively, in the rest frame including the cool core 
(incc hereafter) and correcting for the cool core (cocc here- 
after), respectively. The X-ray luminosity corrected for the cool 
core is computed by assuming a constant value of the surface 
brightness distribution in the cluster core, S'^'^'^iR < Q.lr^^^ ) = 
5x(0.2r^gQ ) for the reason detailed in e.g. Zhang et al. (2007). 

The global properties of R 1504 obey the scaling relations 
of nearby clusters. R1504, for instance, closely matches the 
M -Yx, M - Mgas, M - r, and M - L'=°'='= relations provided 
by e.g. Arnaud et al. (2005, 2007), Vikhlinin et al. (2006), and 
Zhang et al. (2008). This supports the finding that the ICM in 
R1504 is quite relaxed within r^f^\ Regardless of whether the 
high velocity substructure is excluded, R1504 is well within the 
26% scatter of the L'^°^'^ - <t relation for the clusters at similar 
redshifts studied by Zhang et al. (201 la). According to the dy- 
namical mass estimate M2qq'*"'^ and the richness within /"loJ"/'"^, 
R1504 agrees with the mass versus richness scaling relations 
of Johnston et al. (2007) and Hansen et al. (2009), who instead 
used weak-lensing masses. 

5.4. Cluster core 

Strong cool-core systems usually show peaked iron abundances 
in the cluster cores. However, R1504 displays a relatively flat 
radial distribution of iron abundance ranging within [0.2, 0.4] 
times solar abundance (Fig.|6]l. Guo & Mathews (2010) pointed 
out in simulations that AGN outbursts efficiently mix metals in 
the cluster core and may even remove the central abundance 
peak if it is insufficiently broad. 

R1504 is one of nine known clusters hosting radio mini- 
halos. The presence of mini-halos supports the above interpre- 
tation of mixing metals in the cluster core. The BCG is blue 
with a. B - R colour of 0.98, and has strong emission lines. 



The BCG sky position, RA=15*'04"07.573, 5=-02°48'14'.'26 
(J2000), is offset by 4" (14 kpc) from the X-ray flux- weighted 
centroid, which is within the XMM-Newton spacial resolution 
with a 6" full width at half maximum. The BCG harbours a ra- 
dio source with a brightness of 62mJy at 1.4 GHz (Bauer et al. 
2000). Giacintucci et al. (201 1) reveal gas sloshing in the clus- 
ter core, where turbulence may generate particle acceleration 
to form mini-halos. The gas sloshing in the cluster core may 
partially account for the slightly less-concentrated mass distri- 
bution using the X-ray hydrostatic method than that determined 
from the caustic method (Fig.[3]l. 

6. Conclusions 

We have studied one of the most X-ray luminous clusters of 
galaxies in the REFLEX survey, R1504, using XMM-Newton 
imaging spectroscopy, VLT/VIMOS spectroscopy, and WFI 
photometry. 

The mass distribution determined using the hydrostatic 
method based on X-ray imaging spectroscopy agrees with that 
using the caustic method based on optical spectroscopy within 
the Icr uncertainties at most radii, although the former appears 
to be less concentrated than the latter The mass estimate ob- 
tained using the caustic method is 80% of the hydrostatic mass 
estimate at the X-ray mass determined radius r^^-- 

At the overdensity of 500, the dynamical mass estimate is 
72% of the hydrostatic mass estimate. At the overdensity of 
200, the richness-derived mass estimates according to more re- 
cent calibrations of the mass-richness relation (e.g. Johnston et 
al. 2007; Hansen et al. 2009) are about 0.8-1.2 times the hy- 
drostatic mass measurement and 1.0-1.4 times the dynamical 
mass measurement. 

On the basis of the line-of-sight velocities of spectroscopic 
members, our DS test has revealed a relatively high-velocity 
(>1000kms"') group at a projected distance of ~ r^j^- south- 
east of the cluster centroid. The high velocity group was also 
present in the 2-D X-ray maps and weak-lensing shear map. 
The dynamical mass estimate is reduced by 14% when the sub- 
structure of the high velocity group is excluded. This highlights 
the importance of a proper analysis of the structures of evolv- 
ing systems such as cluster of galaxies in order to obtain precise 
mass measurements. 

The cluster R1504 has a rather low stellar-mass fraction 
among nearby clusters (e.g. Biviano & Salucci 2006), although 
there is an agreement given the large scatter for nearby clus- 
ters (e.g. Andreon 2010; Zhang et al. 2011b). Both the intra- 
cluster light and BCG play limited roles in accounting for the 
low stellar-mass fraction. 

Within a projected distance of 700 kpc, the galaxy popula- 
tion is dominated by red galaxies (>70%) based on photometric 
data. Beyond 700 kpc, there is a rapid increase in the fraction 
of blue galaxies. The blue-galaxy fraction reaches about 60% 
at r^Qy ■ and stays almost constant beyond. The dynamical mass 
measurement calculated using blue spectroscopic members is 
1.23 times the value derived using red spectroscopic galaxies. 

Despite the high velocity group in the cluster outskirts and 
the gas sloshing in the cluster core, R1504 appears to be rela- 
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lively relaxed out to /"^go , and obeys the observed scaling rela- 
tions found for similar nearby clusters. 

Acknowledgements. We acknowledge our referee Andrea Biviano, 
who provided insight and expertise that greatly improved the 
work. Y.Y.Z. acknowledges Italo Balestra, Hans Bohringer, Angela 
Bongiomo, Paola Popesso, David Wilman, and Felicia Ziparo for 
useful discussions. Y.Y.Z. acknowledges support by the German 
BMBF through the Verbundforschung under grant 50 OR 1103. The 
XMM-Newton project is an ESA Science Mission with instru- 
ments and contributions directly funded by ESA Member States 
and the USA (NASA). The XMM-Newton project is supported by 
the Bundesministerium fiir Wirtschaft und Technologie/Deutsches 
Zentmm fur Luft- und Raumfahrt (BMWI/DLR, FKZ 50 OX 
0001) and the Max-Planck Society. We have made use of 
VLT/VIMOS observations taken with the ESO Telescope at the 
Paranal Observatory under programme 077.A-0058 and WFI obser- 
vations partially supported by the Deutsche Forschungsgemeinschaft 
(DFG) through Transregional Collaborative Research Centre TRR 33. 
The VLT/VIMOS data presented in this paper were reduced using 
the VIMOS Interactive Pipeline and Graphical Interface (VIPGI) de- 
signed by the VIRMOS Consortium. This research has made use of 
the SIMBAD database, operated at CDS, Strasbourg, France. 

References 

Andreon, S. 2010, MNRAS, 407, 263 

Amal, E. M., Bajaja, E., Larrarte, J. J., et al. 2000, A&AS, 
142, 35 

Amaud, M., Pointecouteau, E., & Pratt, G. W. 2005, A&A, 
441, 893 

Amaud, M., Pointecouteau, E., & Pratt, G. W. 2007, A&A, 
474, L37 

Baade, D., Meisenheimer, K., Iwert, O., et al. 1999, The 
Messenger, 95, 15 

Bajaja, E., Arnal, E. M., Larrarte, J. J., et al. 2005, A&A, 440, 
767 

Balogh, M. L., Navarro, J. R, & Morris, S. L. 2000, ApJ, 540, 
113 

Bartelmann, M. 1996, A&A, 313, 697 

Bauer, F. E., Condon, J. J., Thuan, T. X., & Broderick, J. J. 

2000, ApJS, 129, 547 

Baum, W. A. 1959, PASP, 71, 106 

Becker, R. A., Chambers, J. M., & Wilks, A. R. 1988, Pacific 
Grove, Ca.; Wadsworth Brooks 

Beers, T. C, Flynn, K., & Gebhardt, K. 1990, AJ, 100, 32 
Bell, E. R, & de Jong, R. S. 2001, ApJ, 550, 212 
Bell, E. F, Mcintosh, D. H., Katz, N., & Weinberg, M. D. 
2003, ApJS, 149, 289 

Benin, E., & Arnouts, S. 1996, A&AS, 117, 393 

Biviano, A., Murante, G., Borgani, S., et al. 2006, A&A, 456, 

23 

Biviano, A., & Salucci, P 2006, A&A, 452, 75 

Bohringer, H., Schuecker, P., Guzzo, L., et al. 2001, A&A, 

369, 826 

Bohringer, H., Burwitz, V., Zhang, Y.-Y., Schuecker, P., & 
Nowak, N. 2005, ApJ, 633, 148 

Bohringer, H., Pratt, G. W., Arnaud, M., et al. 2010, A&A, 
514, 32 

Braglia, R, Pierini, D., & Bohringer, H. 2007, A&A, 470, 425 



Butcher, H., & Oemler, A., Jr. 1978, ApJ, 226, 559 
Cappellari, M., & Copin, Y. 2003, MNRAS, 342, 345 
Carlberg, R. G., Yee, H. K. C, Ellingson, E., et al. 1997, ApJ, 
476, L7 

Christlein, D., & Zabludoff, A. I. 2003, ApJ, 591, 764 
Danese, L., de Zotti, G., & di TulHo, G. 1980, A&A, 82, 322 
den Hartog, R., & Katgert, P 1996, MNRAS, 279, 349 
Diaferio, A., & Geller, M. J. 1997, ApJ, 481, 633 
Diaferio, A. 1999, MNRAS, 309, 610 
Diehl, S., & Statler, T. S. 2006, MNRAS, 368, 497 
Dressier, A., & Shectman, S. A. 1988, AJ, 95, 985 
Durret, R, Lagana, T. R, & Haider, M. 201 1, A&A, 529, A38 
Erben, T., Schirmer, M., Dietrich, J. P, et al. 2005, 
Astronomische Nachrichten, 326, 432 

Evrard, A. E., Bialek, J., Busha, M., et al. 2008, ApJ, 672, 122 

Radda, D., Girardi, M., Giuricin, G., Mardirossian, R, & 

Mezzetti, M. 1996, ApJ, 473, 670 

Rujita, Y, & Nagashima, M. 1999, ApJ, 516, 619 

Garim, B., Rumana, M., Rranzetti, P, et al. 2010, PASP, 122, 

827 

Geller, M. J., Diaferio, A., & Kurtz, M. J. 1999, ApJ, 517, L23 
Ghizzardi, S. 2001, in-flight calibration of the PSR for the 
MOSl and MOS2 cameras, EPIC-MCT-TN-011, Internal re- 
port 

Giacintucci, S., Markevitch, M., Brunetti, G., Cassano, R., & 
Venturi, T. 201 1, A&A, 525, LIO 

Gill, S. P D., Knebe, A., & Gibson, B. K. 2005, MNRAS, 
356, 1327 

Girardi, M., Biviano, A., Giuricin, G., Mardirossian, R, & 
Mezzetti, M. 1993, ApJ, 404, 38 

Gladders, M. D., Lopez-Cruz, O., Yee, H. K. C, & Kodama, 
T. 1998, ApJ, 501,571 

Gonzales, A. H., Zaritsky, D., & Zabludoff^, A. I. 2007, ApJ, 
666, 147 

Guo, R, & Mafliews, W. G. 2010, ApJ, 712, 1311 

Hansen, S. M., McKay, T. A., Wechsler, R. H., et al. 2005, 

ApJ, 633, 122 

Hansen, S. M., Sheldon, E. S., Wechsler, R. H., & Koester, 
B.P 2009, ApJ, 699, 1333 

Harrison, E. R., & Noonan, T. W. 1979, ApJ, 232, 18 
Hartmann, D., & Burton, W. B. 1997, Cambridge University 
Press, ISBN 0521471117 

Hildebrandt, H., Erben, T., Dietrich, J. P, et al. 2006, A&A, 
452, 1121 

Hou, A., Parker, L. C, Wilman, D. J., et al. 2012, MNRAS, 
421, 3594 

Johnston, D. E., Sheldon, E. S., Wechsler, R. H., et al. 2007, 

arXiv:0709.1159 

Kaiser, N. 1987, MNRAS, 227, 1 

Kalberla, PM.W., Burton, W.B., Hartmann, D., et al. 2005, 
A&A, 440, 775 

Katgert, P, Mazure, A., Perea, J., et al. 1996, A&A, 310, 8 
Katgert, R, Biviano, A., & Mazure, A. 2004, ApJ, 600, 657 
Kauffmann, G., Colberg, J. M., Diaferio, A., & White, 
S. D. M. 1999, MNRAS, 307, 529 

Lin, Y.-T., Mohr, J. J., & Stanford, S. A. 2004, ApJ, 610, 745 
Mahdavi, A., Hoekstra, H., Babul, A., & Henry, J. P 2008, 
MNRAS, 384, 1567 



10 



Zhang et al.: Probing cluster dynamics in RXCJ1504. 1-0248 via radial and two-dimensional gas and galaxy properties 



Markevitch, M., Forman, W. R., Sarazin, C. L., & Vikhlinin, 
A. 1998, ApJ, 503, 77 

Navarro J. E, Frenk C. S., White S. D. M., 1997, ApJ, 490, 
493 

Odenkirchen, M., Grebel, E. K., Rockosi, C. M., et al. 2001, 
ApJ, 548, LI 65 

Ogrean, G. A., Hatch, N. A., Simionescu, A., et al. 2010, 

MNRAS, 406, 354 

Pisani, A. 1993, MNRAS, 265, 706 

Pisani, A. 1996, MNRAS, 278, 697 

Poggianti, B. M., & Wu, H. 2000, ApJ, 529, 157 

Poole, G. B., Eardal, M. A., Babul, A., et al. 2006, MNRAS, 

373,881 

Popesso, P., Biviano, A., Bohringer, H., Romaniello, M., & 
Voges, W. 2005a, A&A, 433, 431 

Popesso, P., Bohringer, H., Romaniello, M., & Voges, W. 
2005b, A&A, 433, 415 

Popesso, P., Biviano, A., Romaniello, M., & Bohringer, H. 
2007, A&A, 461,411 

Pratt, G. W., Bohringer, H., Croston, J. H., et al. 2007, A&A, 
461,71 

Quilis, v., Moore, B., & Bower, R. 2000, Science, 288, 1617 

Reyes, R., Mandelbaum, R., Hirata, C, Bahcall, N., & Seljak, 

U. 2008, MNRAS, 390, 1157 

Ricker, P M., & Sarazin, C. L. 2001, ApJ, 561, 621 

Rines, K., & Diaferio, A. 2006, ApJ, 132, 1275 

Roettiger, K., Burns, J. O., & Stone, J. M. 1999, ApJ, 518, 

603 

Saintonge, A., Tran, K.-V. H., & Holden, B. P 2008, ApJ, 685, 
L113 

Salpeter, E. E. 1955, ApJ, 121, 161 
Sanders, J. S. 2006, MNRAS, 371, 829 
Schechter, P 1976, ApJ, 203, 297 

Scodeggio, M., Franzetti, P, Garilli, B., et al. 2005, PASP, 
117, 1284 

Snowden, S. L., Mushotzky, R. F, Kuntz, K. D., & Davis, D. 
S. 2008, A&A, 478,615 

Sun, M., Donahue, M., & Voit, G. M. 2007, ApJ, 671, 190 
Urquhart, S. A., Willis, J. P, Hoekstra, H., & Pierre, M. 2010, 
MNRAS, 406, 368 

Verdugo, M., Ziegler, B. L., & Gerken, B. 2008, A&A, 486, 9 
Verdugo, M., Lerchster, M., Bohringer, H., et al. 2012, 
MNRAS, 421, 1949 

Vikhlinin, A., Kravtsov, A., Forman, W., et al. 2006, ApJ, 640, 
691 

Zhang, Y.-Y., Bohringer, H., Finoguenov, A., et al. 2006, 
A&A, 456, 55 

Zhang, Y.-Y, Finoguenov, A., Bohringer, H., et al. 2007, 
A&A, 467, 437 

Zhang, Y.-Y, Finoguenov, A., Bohringer, H., et al. 2008, 
A&A, 482, 451 

Zhang, Y.-Y, Reiprich, T. H., Finoguenov, A., Hudson, D. S., 
& Sarazin, C. L. 2009, ApJ, 699, 1 178 

Zhang, Y.-Y, Okabe, N., Finoguenov, A., et al. 2010, ApJ, 
711, 1033 

Zhang, Y.-Y, Andernach, H., Caretta, C. A., et al. 2011a, 
A&A, 526, A105 



Zhang, Y.-Y, Lagana, T. F, Pierini, D., et al. 2011b, A&A, 
535, A78 

Zibetti, S., White, S. D. M., Schneider, D. P, & Brinkmann, 
J. 2005, MNRAS, 358, 949 

Ziparo, F, Braglia, F G., Pierini, D., et al. 2012, MNRAS, 
420, 2480 



Zhang et al.: Probing cluster dynamics in RXCJ1504. 1-0248 via radial and two-dimensional gas and galaxy properties 1 1 

Table 1. Basic properties of R1504. 



X-ray properties 


JH.E. 
'500 


(1.18 ± 0.03) Mpc 




^"500 


(5.81 ±0.49) X lO^Mo 




rf^ 


(1.76 ±0.05) Mpc 




'"200 


(7.65 ± 0.65) X lO^Mo 




C5OO 


2.93 ±0.14 




Afgas,500(< ''^00') 


(0.907 ± 0.051) X 10'-* Mo 




/gas,500(< 'foo') 


0.156 ±0.015 




Tm.2-0.5)rf^- 


(7.25 ±0.13) keV 




Z(0.2-0.5)/5j|;- 


(0.28 ± 0.03) Zo 




rincc ^H.E.\ 
■^0.I-2.4keV,500'-- '500 
rcocc j.H.E.\ 
-^0.1-2.4keV,500'-- '500 
rincc /-^ r^E.A 
-^bol,500*.- '500 
rcocc .H.E.\ 
-^bol,500^- '500 ' 


(2.13 ± 0.04) X 10-*^^ ergs-' 
(0.61 ±0.04) X 10'*5ergs-' 
(5.48 ± 0.11) X lO-*^ ergs-' 
(1.56±0.11)x 10-*^ ergs-' 




Optical spectroscopic properties with all 53 members 


j-caustic 
'500 


(1.06 ±0.04) Mpc 




v*cau.stic 
'"500 


(4.17 ± 0.42) X 10''*Afo 




caustic 
'200 
yt^caustic 

'"200 


(1.66 ±0.06) Mpc 
(6.43 ± 0.65) X lO'+Mo 




j^causlic/^ H.E.-v 


(4.66 ± 0.47) X lO'^Mo 




C 53m 


(1132±94)kms-' 




^06 
'500,53m 


(1.217 ±0.109) Mpc 




^B06 
'"500.53m 


(6.355 ± 1.71 l)x IO'^'Mg 




Optical spectroscopic properties with 48 members excluding the high velocity group 


0"48m 


(1079 ±97) km S-' 




^06 
500,48m 

'"500.48m 


(1.157 ±0.111) Mpc 
(5.452 ± 1.573) X 10'-* Mq 




Optical spectroscopic properties with all 29 blue members 


0"bluc 


(1170 ± 125)kms-' 




„B06 
'500,bluc 


(1.260± 0.141) Mpc 




'"500,blue 


(7.057 ± 2.373) x IO'^Mq 




Optical spectroscopic properties with all 24 red members 


O-rcd 


(1096 ± 156) kms-' 




„B06 
'500,rcd 


(1.175 ±0.172) Mpc 




MB06 
'"500,red 


(5.722 ± 2.517) X IO'^^Mq 




Photometric properties within X-ray determined r^^g- down to = 22 magnitude 


Galaxy luminosity function 


all members: iJ* = 18.11 ±0.61, a = 
red members: R'^^^ = 18.98 ± 0.34, o-red 


-1.24 ±0.16, CD* = 20.95 ± 12.24 
= -0.86 ±0.19, = 34.87 ± 11.75 


/blue(< ''"oo ) 


0.37 ± 0.06 


A'red(< '■s^oo ) 


180.2 ± 15.5 






283.9 ± 22.4 




iR,bluc(< ''500 ) 


(0.81 ±0.14) X lO'^Lo 




i'R.rcdC^ '500 ) 


(1.77 ± 0.19) X lO'^Lo 




iR,all(< 'i'oo ) 


(2.58 ± 0.24) X lO'^Lo 




M..bluc(< r^^ ) 


(0.31 ±0.07) X lO'-Mo 




M.,red(< rf^ ) 


(3.42 ± 0.66) X lO'-Mo 




M.,all(< rf^ ) 


(3.73 ± 0.66) X lO'^Mo 




,f..500(< r^oo ) 


0.0064 ± 0.0013 




Richness-derived mass measurements down to (R' + 1) magnitude 



(1.851 ±0.082) Mpc 
M^07 (8 94g + 1 190) X lO'-tMc 

Nij!I 92.67 ± 9.63 



r™' (1.651 ±0.077) Mpc 

M™' (6.347 ± 0.887) x 10''' 

83.89 ±9.16 
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E (keV) 



Fig. 1. MOSl (black), MOS2 (red), and pn (green) spectra extracted from the central <0'.5 circle. The lower panel shows devia- 
tions of the models from the data points normalized by the Icr errors. 



> 




r (Mpc) 



R (arcmin) 



Fig. 2. Left: Deprojected temperature measurement versus radius with the parametrized temperature distribution indicated by the 
solid curve and the Icr interval by the dashed curves. The vertical line denotes O.lr^gg . Right: Surface brightness versus projected 
distance from the X-ray flux-weighted centroid with the parametrized surface brightness distribution in solid curve. The upper 
panel shows deviations of the model from the data points normalized by the Icr errors. 
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r (Mpc) 

Fig. 3. Cumulative mass distributions (1) using the X-ray hydrostatic method (filled squares), (2) using the caustic method (black 
solid curve), (3) of the gas mass estimate (black dotted curve), (4) of the stellar mass estimate of the red galaxies (red dash-dotted 
curve), and (5) of the stellar mass estimate of the blue galaxies (blue dash-dotted-dotted-dotted curve). The dashed curves and 
error bars denote lo" intervals. The open star, triangle, diamond, and circle indicate the dynamical mass estimates (M^^) using 
53, 48, 24 red, and 29 blue cluster members according to the Biviano et al. (2006) method. The open square and cross denote the 
richness-derived mass estimates {M'^ and M^^) based on photometric data using the definitions of Johnston et al. (2007) and 
Hansen et al. (2009). 
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Fig. 4. Two-dimensional X-ray spectrally measured temperature (upper left), electron number-density (upper right), entropy 
(lower left), and pressure (lower right) maps using the Sanders (2006) binning scheme. Overlaid circles in blue, red, and green 
denote cluster galaxies that have spectroscopic follow-up data with their clustercentric velocities toward the observer greater 
than 1000 kms"', away from the observer greater than 1000 kms"', and smaller than 1000 kms"', respectively. The radii of 
these circles are proportional to exp[p^g]. The black cross denotes the X-ray flux-weighted centroid and the black dashed circle 
denotes '"^qq'. The white holes mask X-ray detected point sources. The colour scales logarithmically from 4 to 10 in keV for 
temperature, from 0.00005 to 0.12 in cm for electron number density, from 10 to 1800 in keV cm- for entropy, and from 0.0001 
to 0.54 in keV cm"^ for pressure. 
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Fig. 5. Temperature (upper left), electron number-density (upper right), entropy (lower left), and pressure (lower right) distri- 
butions of the spectrally measured 2-D maps using the Sanders (2006) binning scheme are shown as crosses with their local 
regression fits in black curves. The solid and dashed grey curves in the upper left plot show the parametrization of the depro- 
jected temperature profile and its l-cr interval derived in Sect. 12.31 Differential scatter and its 1-cr intervals of the temperature, 
electron number-density, entropy, and pressure fluctuations in the 2-D maps are shown in solid and dashed black curves in their 
upper panels. 
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Fig. 6. Radial distribution of the iron abundance. The vertical line denotes O.lr™ 
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Fig. 7. Left: Histogram of the line-of-sight velocities. Member galaxies of the cluster, member galaxies of the high velocity group, 
and non-member galaxies are in solid (black), dashed (red), and dash-dotted (grey) lines. The arrow shows the cluster redshift. 
Right: Line-of-sight velocity versus projected radius, in which filled blue diamonds and red stars show blue and red member 
galaxies in the caustic (black curve), and open squares show fore- and background galaxies. The grey curve shows the symmetric 
boundary by choosing the minimum of the values on both sides of the caustic boundary relative to the cluster redshift. The cluster 
radii, '"gQQ'*"'^ and '"200^"'^, derived from the mass distribution using the caustic method are shown as solid and dashed vertical lines, 
respectively. 



Zhang et al.: Probing cluster dynamics in RXCJ1504. 1-0248 via radial and two-dimensional gas and galaxy properties 



17 



■ ■ ■ :.-^'n^-^- . . . 



45^00.0 




:3o:o 



1S:O4:0OJ0 



i * 



. •-, I. 




-. *= 




30.0 



I - -1^^ 



;. . s- ::. i fir-- ^■ 




Q3:DOjO 



Fig. 8. WFI fi-band imaging. Overlaid circles in blue, red and green denote cluster galaxies with spectroscopic follow-up data 
that have clustercentric velocities toward the observer greater than 1000 km s', away from the observer greater than 1000 km s^', 
and smaller than 1000 kms ', respectively. The radii of these circles are proportional to expljo^j]. The black cross denotes the 
X-ray flux- weighted centroid and the black dashed circle denotes . 
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Fig. 9. Left: Colour-magnitude diagram with the corresponding red sequence indicated by a solid line and its 3cr — 
0.249 magnitude intervals by dashed lines. The sample is complete down to the curve at ~ 22 magnitude, which varies slightly 
with the B - R colour. We thus limited the analysis to the galaxies brighter than R - 22 magnitude. The dotted line denotes 
the colour cut that is 1.5 magnitude bluer than the red sequence, below which the objects were considered as foreground con- 
tamination. Red stars and blue diamonds indicate red and blue spectroscopic member galaxies regardless of their clustercentric 
distances. The BCG is extremely blue with B - R - 0.98. Right: Cluster galaxy luminosity functions for red (red stars), blue 
(blue diamonds), and all (black circles) members within r^^- derived from the X-ray hydrostatic mass estimate. The membership 
is statistically determined using control fields with the background contamination statistically subtracted. The best-fit Schechter 
(1976) functions are shown in red and black curves for the red and all members, respectively. The inset shows the 68% confidence 
region for the R* and a parameters of the best-fit Schechter (1976) functions. No fit was possible for the blue members. 




Fig. 10. Left: Number density profiles with their best-fit NFW models for red member galaxies (red triangles, solid curve) and all 
members (black circles, dashed curve), respectively. We obtained reduced j^'^ = 1.03, 2.24 for the red and all member galaxies, 
which indicates that the red galaxies follow well the theoretical profile and the blue galaxies may account for the deviation of 
the fit for all member galaxies. The bin size is 300 kpc. Right: Ratio of blue to all member galaxies as a function of projection 
distance from the X-ray flux-weighted centroid. The vertical line remarks r^,^'- 
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Appendix A: X-ray 2-D spectrally measured maps 
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Fig. A.l. Two-dimensional X-ray spectrally measured temperature (upper left), electron number-density (upper right), entropy 
(lower left), and pressure (lower right) maps using the Cappellari & Copin (2003) binning scheme. Overlaid circles in blue, red, 
and green denote cluster galaxies that have spectroscopic follow-up data with their clustercentric velocities toward the observer 
greater than 1000km s"', away from the observer greater than 1000km s ', and smaller than 1000km s ', respectively. The radii 
of these circles are proportional to exp[ppg]. The black cross denotes the X-ray flux-weighted centroid and the black dashed 
circle denotes /"^qo ■ white holes mask X-ray detected point sources. The colour scales logarithmically from 4 to 10 in keV 
for temperature, from 0.00005 to 0.12 in cm""* for electron number density, from 10 to 1800 in keV cm' for entropy, and from 
0.0001 to 0.54 in keV cm"^ for pressure. 



